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Abstract 

We have observed a bimodal behaviour of the distribution of the asymmetry 
between the charges of the two heaviest products resulting from the decay of the 
quasi-projectile released in binary Xe+Sn and Au+Au collisions from 60 to 100 
MeV/u. Event sorting has been achieved through the transverse energy of light 
charged particles emitted on the quasi-target side, thus avoiding artificial correla- 
tions between the bimodality signal and the sorting variable. Bimodality is observed 
for intermediate impact parameters for which the quasi-projectile is identified. A 
simulation shows that the deexcitation step rather than the geometry of the colli- 
sion appears responsible for the bimodal behaviour. The influence of mid-rapidity 
emission has been verified. The two bumps of the bimodal distribution correspond 
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to different excitation energies and similar temperatures. It is also shown that it 
is possible to correlate the bimodality signal with a change in the distribution of 
the heaviest fragment charge and a peak in potential energy fluctuations. All to- 
gether, this set of data is coherent with what would be expected in a finite system 
if the corresponding system in the thermodynamic limit exhibits a first order phase 
transition. 



1 Introduction 

Many experimental features suggestive of the liquid-gas phase transition of nu- 
clear matter have been observed in the Fermi energy regime of nucleus-nucleus 
collisions. Among these features are: abnormal partial energy fluctuations [1,2], 
charge correlations [3], a double peaked distribution of an order parameter 
("bimodality") [4,5], fluctuation properties of the heaviest fragment size [6, 7], 
Fisher scaling[8], vaporisation [9], flattening of the caloric curve[10,ll]. In the 
present paper, we concentrate on the bimodality signal. 

The discontinuity of the order parameter at a first order phase transition is 
expected to be replaced in a finite system by a bimodal distribution of the 
order parameter close to the transition point [12]. If the order parameter is one 
dimensional [15] and if the transition has a finite latent heat, the bimodality 
should be observed in the canonical ensemble, the transition temperature being 
denned as that at which the two peaks have the same height [13]. In the 
fragmentation transition case, the size of the heaviest cluster produced in 
each event is an order parameter in many different models [14] including the 
lattice-gas model of the liquid-gas phase transition[15]. In this picture, when 
the size of the heaviest fragment is large, the system is mostly on the liquid-like 
side, whereas if it is small one is mostly dealing with a gas-like behaviour. 

An example is given in figure 1 obtained in the lattice-gas approach (see 
for instance ref[16] [13]). For temperatures outside of a small range around 
the transition temperature, the size of the largest fragment exhibits single- 
humped distributions whereas two peaks of the same height are obtained at 
the transition temperature value. This result which is obtained for an equili- 
brated system in the canonical isobar ensemble remains approximately valid 
even if part of the available energy is not thermalized. This is shown in figure 

2 in which an aligned momentum has been randomly added for each parti- 
cle, which simulates the dynamical effects corresponding to the memory of 
the aligned beam momentum in the entrance channel. Several percentages 
of aligned momenta have been considered ranging from zero percent (purely 
thermal situation) to 100 percent (same amount of aligned and thermalized 
momentum). The bimodality signal is robust in the sense that it is observed 
even if a sizeable amount of the total available momentum is still aligned along 
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the beam direction, with only a slight change of the deduced transition tem- 
perature. Only if the aligned momentum is as important as the thermal one 
(lower right in the figure) the distribution is weakly bimodal and mimics a flat 
distribution expected for a continuous transition [12, 14]. 

Guided by these considerations, preliminary evidences of bimodal behaviours 
in some static observables have been obtained in previous data concerning the 
systems Ni+Au[18], Ni+Ni[19], and Xe+Sn[20]. In none of these cases, a clear 
bimodality was observed in the distribution of the heaviest fragment charge. A 
possible reason is that references [18,19,20] concern central collisions. In such 
cases, the event selection is closer to a microcanonical sorting, with a selection 
of a narrow range of excitation energy: if the transition is characterized by a 
finite latent heat, the two peaks that can be interpreted as precursors of the 
two coexisting phases [16] should be associated to two different energies. The 
signal which would naturally appear if the system was in contact with a heat 
bath is expected to be quenched if the experimental sorting constrains too 
strongly the deposited excitation energy. 
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Fig. 1. Size distribution of the largest 
product obtained with the lattice-gas 
model for an equilibrated system of 216 
particles in the canonical isobar ensem- 
ble at subcritical pressure[13]. e repre- 
sents the lattice coupling. The three 
plots correspond to temperatures re- 
spectively below, at and above the 
transition one. 




Fig. 2. Similar to figure 1 but an aligned 
momentum component has been added 
for each particle. The added aligned 
momentum is 0, 10, 50 or 100 percent 
of the thermalized one. The last insert 
corresponds to the 100 percent case 
at three different temperatures. Such 
calculations which include an out-of 
equilibrium component are presented 
in reference [13]. 



In order to mimic a canonical sorting in the data selection, one can select semi- 
peripheral collisions which are mostly binary. In these collisions, it is possible 
to isolate a quasi-projectile (QP) and a quasi-target (QT). The fluctuations 
in the sharing of the dissipated energy between these two partners may be 
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sufficient to explore the two phases. In the present work, we focus on such 
peripheral and semi-peripheral collisions. 



2 The experiment 

The data have been obtained by the Indra-Aladin collaboration (Indra at 
GSI). The Indra set up is described in references [2 1,2 2], while details of the 
measurements performed at the GSI laboratory and of the analysis and cal- 
ibration procedures may be found in [23,24,25,26]. The studied systems are 
Xe+Sn and Au+Au from 60 to 100 MeV/u. The events were registered if at 
least four modules fired. In this section the question of the event sorting is 
first addressed. Then we will check to which extent it is possible to isolate 
properly the QP from the QT and any mid-rapidity emission. 

2.1 Event sorting 

In order to separate the QP and QT contributions, the momentum tensor 
calculated from Z > 3 products has been diagonalized event by event. The 
main axis of the tensor is in most cases close to the beam axis. The particles 
and fragments attributed to the QT (resp. QP) are those which are backward 
(resp. forward) emitted along the main tensor axis. Event sorting has been 
achieved from the transverse energy Etrans of light charged particles (LCP) 
on the QT side calculated in this frame. A selected Etrans value corresponds 
to a selected violence of the collision. 

The Etrans sorting has the advantage of minimizing autocorrelations between 
the sorting variable and observables, since the sorting particles are only loosely 
correlated to the particles and fragments considered to get the QP properties. 
From the experimental point of view, the detection efficiency is good for LCP 
in the whole space since the LCP detection thresholds are small. Hence, the 
LCP transverse energy on the QT side was well measured. On the other hand, 
the detection efficiency on the QP side (velocities exceeding the cm. velocity) 
is good for both the LCP and the fragments since all the laboratory veloci- 
ties are far above the detection thresholds for any product. In order to ensure 
an almost complete QP reconstruction, any event for which the forward de- 
tected charge was below 80 percent of the projectile charge has been rejected. 
This selection rejected mainly peripheral events for which the QP residue was 
emitted in the forward detector opening. 

The Etrans distributions, normalized to the size of the system and to the 
laboratory bombarding energy in MeV/u, for different systems are shown in 
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normalized Etrans 

Fig. 3. Etrans distributions for the two studied systems at the three indicated 
energies. For the abscissa Etrans has been divided by the total mass number of the 
system (projectile + target), and by the incident energy in MeV/u. 

figure 3. They are similar and they extend on the same abscissa range whatever 
the system and the incident energy. The differences between the systems are 
mainly due to the loss of efficiency for very peripheral collisions when the 
quasi-projectile residue escapes in the detector forward hole (these events are 
rejected because of the completeness detection criterion). This efficiency loss is 
stronger for the lighter system because QP fission has a negligible probability 
in the Xe mass range: the QP loss in the forward hole is indeed reduced 
when fission occurs because the fission fragment emission angles are generally 
pushed away from very forward directions. 

The similarity of the various normalized- -Etrans distributions indicates that 
the energy damping (or the stopping power) does not depend strongly on the 
mass number nor on the available energy in this bombarding energy region. 
It is in agreement with the conclusions of reference [2 7]. The largest values 
reached by the normalized Etrans are about 0.04. They can be compared 
with the limiting value of 1/12 ~ 0.08 which would correspond (for such 
symmetrical systems) to a complete fusion followed by vaporization (with 
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Fig. 4. Alpha particle centre-of-mass v. 



par 



-'per 



vpar (cm/ ns) 

plots for nine Etrans bins covering 
the whole Etrans distribution of figure 3. Xe+Sn system at 80 MeV/u.The width 
of each Etrans bin is 100 MeV i.e. 0.4 MeV/u if Etrans is normalized on the total 
system mass as in figures 11 and 12 

only LCPs as decay products). 



2.2 Kinematical separation between QP and QT decay products 



This paper is devoted to the QP decay study. Hence, it must be verified to 
which extent the QP and QT decay products can really be separated. This 
question can be answered looking at figure 4 in which typical v par — v perp 
invariant cross section plots are shown for alpha particles emitted in 80 MeV/u 
Xe on Sn collisions. v par and v per are the parallel and perpendicular velocity 
components in the reaction centre of mass frame. Similar distributions are 
obtained for the other systems considered in this paper. The events have been 
sorted in nine Etrans bins. The distributions clearly exhibit two sources for 
peripheral and semi-peripheral collisions at least up to the fourth bin. For the 
highest Etrans values the distinction between a QP and a QT source becomes 
somewhat artificial. 

A careful observation of figure 4 indicates that alpha particles are not isotrop- 
ically emitted from the QP and QT. Instead more alphas are emitted at mid- 
rapidity in the neck region. This result is a general feature in intermediate 
energy nucleus-nucleus collisions[28,29,30] and one may worry about the in- 
fluence of a mid-rapidity emission on the results discussed in this paper. This 
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will be discussed in section 5. 
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Fig. 5. Correlation between the second 
largest charge Zmax2 and the largest 
charge Zmax\ for the Au+Au system 
at 80 MeV/u. 
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Fig. 6. Correlation between the second 
largest charge Zmax2 and the largest 
charge Zmax\ for the Xe+Sn system 
at 80 MeV/u. 



3 QP observables 

Following the discussion of the introduction, a natural choice of the order 
parameter is the size (charge) of the largest fragment Zmax\. The asymme- 
try Z asym between the two heaviest products Zmax\ and Zmax2 (see below 
equation 1 for its precise definition) will also be studied. 

In figures 5 and 6, the correlations between Zmax\ and Zmax2 are shown for 
the Au+Au and Xe+Sn systems, respectively, at 80 MeV/u. All the retained 
(complete) events are included in the picture. For the Au+Au system, the 
binary fission of the quasi-gold nucleus can be recognized for Zmax\ and 
Zmax2 values around 40 charge units. This fission contribution is negligible 
for the Xe case. For both reactions one recognizes also two dominant classes of 
events. Residue like events correspond to a large Zmax\ and a small Zmax2- 
Multifragmentation events correspond to small values for both Zmax\ and 
Zmax2- An important feature is that there is not a continuous evolution from 
one class of events to the other. Instead, they are widely separated showing 
that they correspond to two clearly distinct mechanisms. 

For Au QP, the presence of the fission channel perturbs the evolution of Z asym 
since this variable exhibits similar values for multifragmentation and binary 
fission. In order to avoid this difficulty, we applied a specific treatement to 
the fission events observed for the Au case. In effect, they belong also to the 
normal compound nucleus decay event class and can be classified together 
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with the residue events. For this purpose, we defined Z max as the sum of the 
two fission fragment charges when fission was recognized. Fission events have 
been selected from the charge product of the two fission fragments: charge 
product exceeding 900 [1] (see also figure 5). For these events, Z max is the sum 
of the two heaviest products and we defined Z sec as the third heaviest product 
charge. For the other Au events and for any event in the Xe case (for which the 
fission contribution is negligible), one has Z max = Zmax\ and Z sec = Zmaxi 
. Then the asymmetry parameter is defined by: 
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In order to avoid spurious peaks due to integer number effects we have added 
a random number between -O.i 
exceed slightly the range [0-1]. 



a random number between -0.5 and +0.5 to Z max and Z sec . Thus Z asym can 
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Fig. 7. Z n 



J asym 



plots for the 



nine Etrans bins. Au+Au system at 80 
MeV/u. 



Fie. 8. Z n 



J asym 



plots for the 



nine Etrans bins. Xe+Sn system at 80 
MeV/u. 



4 Bimodality observation 



Let us consider now the evolution of the correlation between Z max and Z asym 
when the events are sorted according to the transverse energy Etrans defined 
in section 2.1. Figures 7 and 8 correspond to the systems Au+Au and Xe+Sn, 
respectively, at 80 MeV/u. They show the correlations between the asymmetry 
variable Z asym and the heaviest product charge Z max for the nine Etrans bins 
already used in figure 4. Two behaviours are clearly evidenced in the figures. In 
the first bins, the heaviest decay product is generally a residue (large Z max and 
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Fig. 9. Zmax distributions for the 

9 Etrans bins. Xe+Sn system at 80 
MeV/u. 
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Fig. 10. Z asym distributions for the 

9 Etrans bins. Xe+Sn system at 80 
MeV/u. 



Z aS ym value close to unity). By contrast, the most violent (central) collisions 
correspond to small values of Z max and Z asym associated with multi-fragment 
emission (such observations are quite in agreement with the results of ref [31]). 
In between, in the third bin, both behaviours are observed: the distribution is 
bimodal. Similar results have been obtained by Ma et al[32] on a lighter system. 
Here we have such an observation on medium and large mass systems. With 
the definition of Z, max as the sum of the two fission fragment charges (see the 
previous section), the fission event behaviour is quite similar to the residue 
one. We have checked that, if one does not reconstruct the parent nucleus 
charge before fission, these events correspond to an additional bump in figure 
7 widely separated from the residue and the fragmentation bumps. It does not 
modify the bimodality behaviour observed in the third bin. 

Although the event probability distribution in the third bin exhibits a mini- 
mum in the representation of figures 7 and 8, this minimum (convex shape) is 
not apparent on the projections on the two shown in figures 9 and 10 

for the Xe+Sn system at 80 MeV/u. Nevertheless, these figures indicate that 
the most probable values of the Z max or Z asym distributions change dramati- 
cally between the third and the fourth bins. This evolution is summarized in 
figures 11 and 12 for the two studied systems and the three bombarding ener- 
gies. In these figures, the abscissa is the transverse energy sorting parameter 
Etrans divided by the system size (sum of the projectile and target masses). 

The abrupt change of the order parameter together with the wide and al- 
most structureless distribution at the transition point shown in figures 7-10 
are rather suggestive of a continuous (second order) transition point. However, 
as we have discussed in section 1, such a behaviour is also consistent with a 
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first order transition polluted by strong out of equilibrium effects, or impor- 
tant deformations due to sorting constraints [15]: if the sorting variable is too 
strongly connected to the order parameter, its bimodal character is naturally 
suppressed. 
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Fig. 11. System Au + Au at 60, 80, 
100 MeV/u (respectively open points, 
stars and black points). Evolution with 
Etrans of the most probable value of 
the Z asym distributions. In this figure, 
Etrans has been divided by the total 
mass number of the system (projectile 
+ target). It has not been divided by 
the incident energy as in figure 3. 
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Fig. 12. Similar to figure 11 but for the 
systems Xe + Sn at 65, 80 and 100 
MeV/u. The abscissa corresponding to 
the transitions are similar for Au + Au 
and Xe + Sn system. This scaling in- 
dicates that the transverse energy per 
nucleon is a key parameter. 



5 Is bimodality due to dynamics? 



In comparing figures 11 and 12, one may conclude that the transverse energy 
per nucleon at which the transition occurs does not depend on the system 
size. A similar kind of universality has also been pointed out in reference [33]. 
Such a result is expected if the transition reflects the energy density involved 
in the reaction. The transition energy for a definite system evolves when the 
bombarding energy increases: it is roughly proportional to the incident energy. 
This result is not expected if Etrans reflects the thermal energy stored in the 
quasi-projectile. 

However, Etrans reflects both the dissipated energy and preequilibrium ef- 
fects. When one increases the bombarding energy, the preequilibium contri- 
bution increases thus inducing a larger value of Etrans to reach the same 
thermal energy. A similar behaviour has been evidenced from the calculations 
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Fig. 13. Angular distributions of the heaviest QP fragment in the QP frame. The 
abscissa is the cosine of the angle between the QP direction in the centre of mass 
system and the Zmaxl fragment direction in the QP frame. The events have been 
sorted according to the Etrans variable. Most of the heaviest fragments are forward 
emitted. Xe+Sn system at 80 MeV/u. 

discussed in figure 2. To test such an interpretation, we have to introduce an 
additional selection of the events in order to decrease the influence of preequi- 
librium effects. Section 5.3 is devoted to the corresponding data analysis. 

One may also consider that the transition (where bimodality is observed) cor- 
responds to an intermediate impact parameter which would be about the same 
whatever the bombarding energy is. If for this centrality a massive fragment 
emission from the neck region sets in, the sudden change of the Z max distri- 
bution could be related to a sudden decrease of the QP source mass, i.e. to a 
geometrical effect. This interpretation is discussed in sections 5.1 and 5.2. . 

5. 1 The role of midrapidity emission. 

Dynamical effects were clearly recognized in nucleus-nucleus collisions in the 
intermediate energy regime. We already noticed in figure 4 the possible role 
of neck emission which has been extensively studied in reference [24,34,35]. A 
generalization of the results of reference[34] is shown in figure 13 where the 
angular distributions of the heaviest fragment (Zmaxl) are drawn for various 
Etrans bins. The abscissa is the cosine of its emission angle in the QP frame. 
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For this analysis, the QP centre of mass has been calculated including all 
products with Z > 2 emitted forward in the ellipsoid frame in order to be able 
to treat events with a residue and no intermediate mass fragment (IMF: Z > 
3). The distributions would be flat for isotropic emission which is not the case. 
Mostly, the heaviest fragment is forward emitted in the QP frame. This feature 
is a clear entrance channel memory: full equilibrium has not been achieved 
at least for most events and light IMF or LCP are emitted at mid-rapidity 
between the QP and QT. What is the role of this mid-rapidity emission? Does 
it contribute significantly to the sudden increase of the fragment multiplicity 
when passing from the residue-like class to the multifragmentation one? 

Charge, Z=1, Z = 2, Z = 3-5, Z = 6-9, Z = 9-14, Z ge 1 5, Vsym ge .8 
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Fig. 14. Charge density distributions for various IMF charges as a function of the 
parallel velocity to the beam in the cm. frame in c units. The full and open circles are 
the distributions for Etrans bins 3 and 9 respectively. They have been normalized 
for zero velocity. The black squares correspond to the difference between the two 
previous distributions. For the third bin the events with Z asym larger than 0.8 have 
been selected. System Au + Au at 100 MeV/u. 

In order to evaluate the role of mid-rapidity we performed an analysis based 
on the hypothesis that mid-rapidity emission comes mainly from the overlap 
region between the projectile and the target and that its kinematical properties 
are the same for all the impact parameters. This assumption is in agreement 
with the basic hypothesis of the HIPSE model [37] which is discussed in section 
5.2. 



Central collisions are those for which the overlap between the projectile and 
the target is the largest: the overlap or participant zone is in this case the 
unique source for particle emission. The corresponding velocity properties have 
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been assumed to be the reference properties of mid-rapidity emission at any 
impact parameter (similar parallel velocity distribution shapes for any impact 
parameter). They have been summarized in plotting the charge density p as a 
function of the velocity parallel to the beam v par in the centre of mass frame. 

The Au + Au system at 100 MeV/u will be taken as an example. Central 
collisions have been associated with the largest Etrans bin. As a matter of fact, 
one may see in figure 4 that the vpar—vper plots for the ninth bin is dominated 
by mid-rapidity emission with a negligible contribution of QP or QT emissions. 
Now, we focused on the third Etrans bin where the bimodality signal is seen. 
In figures 14 (resp. 15) we have plotted the p(vpar) distributions for several 
IMF charges and for Z asym larger than 0.8 for residue-like selection (resp. lower 
than 0.2 for multifragmentation). The black dots correspond to the data. The 
open ones are the corresponding reference distribution obtained for central 
events (ninth bin). The two distributions have been normalized at zero velocity 
which means that we assumed that the contribution of the QP at mid-rapidity 
was negligible. The difference between the two distributions (black squares) is 
attributed to the sequential QP decay. In figure 14 the statistics is very poor 
for Z between 6 and 14 because very few IMF are emitted when a residue 
is released. On the contrary, a larger number of IMF is obtained when one 
selects multifragmentation (figure 15). The resulting QP contributions (black 
squares) for all charges are symmetric around the QP velocity derived from 
the bottom right frame as it is expected for late sequential decay. This means 
that any QP decay contribution with an influence from the QT partner, which 
would show a higher yield of particles in the backward direction [38, 39], has 
been included in the mid-rapidity contribution. The contribution attributed 
to the QP in this analysis is hence underestimated. 
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10.4 


< 0,2 


Total 


19.95 


5.37 


14.57 


68.3 


41.1 


27.2 



Table 1 

Repartition between the equilibrated QP and the mid-rapidity of the total available 
charge forward emitted in the centre of mass. The first three lines are for the residue- 
like events. The three last ones correspond to multifragmentation events. They all 
belong to the third Etrans bin for which bimodality is observed. M m j_ mp , Mqp, 
Mtot are the multiplicities respectively associated to mid-rapidity, QP emission, and 
any emission for LCP or IMF or any charged product. Z m j_ rap , Zqp, Z tot are the 
corresponding charges. Au + Au system at 100 MeV/u. 
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Fig. 15. Similar to figure 14 but the events selected for the third bin are multifrag- 
mentation events (Z asym lower than 0.2). 

One may now deduce from the above treatment the minimum QP charge 
or the maximum mid-rapidity contribution in the third Etrans bin of inter- 
est, i.e. when bimodality is observed. The results are indicated in table 1. 
The mid-rapidity contribution is never dominant and the mean IMF mul- 
tiplicity attributed to the QP increases from 1.36 to 3 from residue like to 
multifragment events. This result seems to indicate that the bimodality signal 
(disappearance of the residue and strong increase of IMF production) is not 
dominated by mid-rapidity emission. 



5.2 Comparison with a dynamical model. 



The role of dynamics has been tested with the help of the event generator 
HIPSE[37] in which the geometry of the collision is explicitly taken into ac- 
count. For any impact parameter the overlap region leads to mid-rapidity 
emission, and contributes to QP excitation from nucleon-nucleon collisions. 
This is thus particularly suited to explore the puzzling property observed in 
figures 11 and 12: the bimodality transition is always observed for similar 
impact parameters corresponding to a Etrans value which increases roughly 
proportionally with the beam energy. The parameters of the model have been 
adjusted in order to reproduce most of the experimental inclusive distribu- 
tions (charge, angular and energy distributions) for Xe+Sn reactions between 
25 and 50 MeV/u. The global experimental features are well reproduced by 
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the model, as shown in figure 16 for the Xe + Sn system at 80 MeV/u. This 
figure is similar to the experimental one (figure 8) and one recognizes bimodal 
distributions in the third and the fourth bins in very good agreement with 
the data. In the model, it is possible to go further and search the reason why 
the residue production vanishes when Etrans is increased. The reason for this 
change is not directly the fact that a decrease of the impact parameter leads 
to a smaller QP spectator. This conclusion may be drawn from figure 17 in 
which is shown the correlation between the initial hot QP size and Z asym for 
the nine Etrans bins. It turns out that the initial QP size is always large. 
This is especially the case for bins 3 and 4, i.e. in the region where bimodality 
is observed in figure 16. This result means that it is the decay step which is 
responsible for the disappearance of the residue. It is found from the model 
that both the excitation energy and the angular momentum of the QP be- 
come quite large (3-5 MeV/u and 70h) when the abrupt change in the decay 
mechanism is observed. In reference [40], by using a sequential decay code, 
the authors conclude that the angular momentum is the key parameter. To 
progress on the understanding of the nature of the transition it will be nec- 
essary to perform similar simulations with a non- sequential decay code and, 
from the experimental point of view, to compare systems with very different 
values of angular momentum but similar excitation energies in the range of 
3-5 MeV/u. 
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Fig. 16. Z asym — Z max plots for the 9 Etrans bins as it is calculated by the HIPSE 
event generator. Xe+Sn system at 80 MeV/u. 
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Fig. 17. Z asym — Zqp plots for the 9 Etrans bins as it is calculated by the HIPSE 
event generator. Xe+Sn system at 80 MeV/u. 
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Fig. 18. Zmax distributions for the 9 
Etrans bins when the "less dynamical 
events" (see text) are selected. Xe+Sn 
system at 80 MeV/u. 
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Fig. 19. Zasym distributions for the 9 
Etrans bins, when the "less dynamical 
events" (see text) are selected. Xe+Sn 
system at 80 MeV/u. 



5.3 Test of preequilibrium effects 



As shown in section 5.1, the angular distributions of the heaviest fragment, 
(Zmaxl), are forward peaked (figure 13), which indicates a memory of the en- 
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trance channel. Nevertheless, for some events, the heaviest fragment is emitted 
backward in the QP frame. The memory of the entrance channel is weaker for 
these events which will thus be labelled "less dynamical events". In figures 
18 and 19, they have been isolated in selecting cos(9) values lower than —0.4 
(see figure 13). By comparison with those obtained when this condition is not 
required (figures 9 and 10), two facts may be emphasized: firstly the bimodal 
behaviour is better evidenced when the "less dynamical events" are selected. 
This is especially the case for the Z asym variable which exhibits clearly a dou- 
ble humped distribution in the third Etrans bin. On the other hand, the 
transition takes place for a lower Etrans value for these selected events. All 
these features are expected if bimodality reflects a thermal behaviour. Indeed, 
lattice-gas calculations (figure 2) indicate that the bimodality signal is attenu- 
ated by dynamical effects; the observed difference between the two samples of 
events may be explained in this framework. Moreover, when dynamical effects 
are present, we can expect the value of Etrans at the transition point to be 
increased since the contribution of a non-relaxed energy component does not 
modify significantly the bimodality signal. 



6 Transition excitation energies and temperature 

To progress in the understanding of the bimodality signal, in this section, we 
try to connect its properties to signals of liquid-gas transition. If fragmentation 
can be associated to such a transition, then a bimodality should be observed 
in the probability distribution of the system density which is the order pa- 
rameter of the transition. Two measurable variables strongly correlated to the 
density are the size of the heaviest fragment and the energy: both should have 
bimodal probability distributions in the liquid-gas transition as illustrated in 
the framework of the lattice-gas model in figure 20. On the experimental side, 
the matter density is not accessible but the excitation energy can be obtained 
from calorimetry as described in section 6.1. 

6. 1 QP reconstruction 

In order to be able to characterize the QP, one has to identify its centre 
of mass and its excitation energy. As discussed in section 2.1, any product 
emitted forward in the ellipsoid frame has been attributed to the QP. In order 
to minimize the contribution of midrapidity, we have calculated the origin 
of the QP frame from the IMF only (Z > 3) at variance from section 5.1. 
Then, we have estimated the deposited energy from calorimetry in including 
all the IMF and only the LCP forward emitted in the QP frame. The LCP 
contribution has been doubled to account for backward emitted LCP. This 
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Fig. 20. Correlation between the largest 
fragment size and the internal exci- 
tation energy from lattice-gas calcula- 
tions close to the transition tempera- 
ture (see ref [13]). The projections are 
also shown. Bimodality is observed for 
both variables. 




Fig. 21. Excitation energy (MeV/u) 
distributions for various Etrans bins. 
The events have been separated in two 
classes corresponding to the 2 bimodal- 
ity solutions: Z asym > 0.7 (full line 
histograms) and Z aS y m < 0.3 (dashed 
line histograms). Xe+Sn system at 80 
MeV/u. 



procedure permits to avoid to a large extent LCPs resulting from mid-rapidity 
emission. The charge of the original QP has been determined by adding the 
detected product charges. Its mass has been taken as that giving the same 
N/Z ratio as the projectile. The masses of the fragments (for Z exceeding 4) 
have been given by the EPAX formula[42] and the neutron number has then 
been obtained from the difference between the initial QP mass and the total 
masses attributed to the detected products. Then the excitation energy has 
been calculated from: 



MlMF M L CP 

E*= E T *+ E Tj + M neut {T n ) - Q 
i=i j=i 



(2) 



where the values are the kinetic energies for each product, Mjmf, Mlcp 
and M neut are the multiplicities of IMF, LCP and neutrons and Q the mass 
balance. The average kinetic energy of neutrons (T n ) has been obtained by 
means of an effective temperature [45,46]. 



6.2 Excitation energy distributions 



In figure 21, the excitation energy distributions are shown for various Etrans 
bins and for the Xe + Sn system at 80 MeV/u. For each Etrans bin, we 
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have separated the events corresponding to Z asym larger (resp. lower) than 
0.7 (resp. 0.3). The two distributions have similar amplitudes in the third 
bin in which the bimodality signal is observed (see figure 8) but they do 
not coincide. Even if excitation energies as high as 15 MeV/u are certainly 
not realistic, we can see that a larger mean excitation energy is measured 
for the multifragmentation case (low Z asym ) in qualitative agreement with the 
calculation results of figure 20. Many reasons can be invoked to explain the fact 
that these two distributions overlap more strongly than the model calculation 
results: the sorting is not a canonical one; moreover, the calorimetry method 
is far from being perfect due to the doubling procedure and to the estimated 
neutron correction [4 7]. 
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Fig. 22. Temperatures measured as a function of Etrans for residue and mul- 
tifragment events. Xe+Sn system at 80 MeV/u. The black (resp. open) points 
correspond to two thermometers: slope parameter of protons (resp. double iso- 
topic ratio method). The error bars are statistical. The transition region is around 
Etrans/A = lMeV/u (figure 11). 

6. 3 Temperatures 

One may now try to measure the involved temperatures. If bimodality has 
a thermal origin, one would expect to observe similar temperatures for both 
classes of events: residue-like events and multifragmentation. 

Temperature measurement is a difficult task because the various possible ther- 
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mometers are not equally suited to the different energy regimes. The kinetic 
energy slope method is expected to be adapted when a residue is released, but 
not in the multifragmentation case because the particle source is not unique. 
The double isotopic ratio is better suited in this second case because the sys- 
tem can explore a small density region [48] and because the IMF multiplicity 
is significant. In figure 22, we indicate the corresponding results for the four 
first bins of figure 21. The abscissa is the transverse energy normalized to the 
system size and the ordinate indicates the measured temperatures. The dou- 
ble isotopic ratio method (He-Li nuclei) has been used for the events selected 
with Z asym lower than 0.3 and the corresponding results have been corrected 
for side feeding effects[48]; the slope of proton spectra (forward emitted in 
the QP frame) has been used for the residue events (Z asym larger than 0.7). 
The points with low statistics have not been plotted. In the transition energy 
region (0.8 — 1.2 MeV/u), both groups exhibit similar temperature values as 
expected if bimodality has a thermal origin. 

The events belonging to the two regions of the bimodal distribution are thus 
found to correspond to different excitation energies but to similar tempera- 
tures. It is what is expected for a first order phase transition for which the 
energy jump between the two solutions reflects the latent heat. This result 
is hence coherent with a thermal origin of multifragmentation but cannot be 
considered as a strong evidence because the compared temperatures have been 
obtained with two different thermometers. 



7 Bimodality and other possible phase transition signals 

As mentioned in the introduction, several signatures of phase transition have 
been proposed and experimentally observed. Due to the various experimental 
conditions, each of these signals present some weaknesses. Thus the occur- 
rence of a phase transition would be strongly reinforced by a simultaneous 
observation of several signatures. Such is the aim of this section. 

7.1 Relation between bimodality and A-scaling 

A-scaling is a signal proposed in the framework of the universal fluctuation 
theory [14]. It may be used to distinguish between different phases and to 
identify critical points. A-scaling is observed for an order parameter m if the 
probability P/v( m ) for observing a given value of m can be written with the 
following formula 

< m > A P N (m) = f(z A ) = f( y — - r l) (3) 
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Fig. 23. Normalized Z max distributions 
according to equation 3 for the sys- 
tem Xe+Sn at 80 MeV/u with A=0.5. 
Various symbols correspond to various 
Etrans bins. 
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Fig. 25. Delta scaling for the sys- 
tem Xe+Sn at 80 MeV/u with A=0.5. 
Only residue like events (Z asym larger 
than 0.7) have been retained. Various 
symbols correspond to various Etrans 
bins. Only the Etrans bins for which 
the residue contribution is dominant 
have been retained. 



Fig. 24. Similar to figure 23 with A=l. 
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Fig. 26. Delta scaling for the sys- 
tem Xe+Sn at 80 MeV/u with 
A=l. Only multifragmentation events 
(Z aS ym lower than 0.3) have been re- 
tained. Various symbols correspond to 
various Etrans bins. Only the Etrans 
bins for which the multifragmentation 
contribution is dominant have been re- 
tained. 



in which /(-2a) is a universal function which does not depend on the size N 
of the system. For a system exhibiting two phases, the A value is expected 
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to be close to 0.5 for an ordered phase and to 1 for a disordered phase. The 
transition point between the two phases is first order if, at this point, the 
distribution of m is bimodal, and second order if it obeys A=l scaling with 
a non-Gaussian tail [14]. Several data have been analyzed in this context [6,7]. 
In reference [7] A-scaling was evidenced for central collisions for the Xe + Sn 
and Au + Au systems of interest in the present paper. Those results indicate 
that the size of the largest fragment is an order parameter. The ordered (resp. 
disordered) behaviour is observed below (resp. above) an incident energy of 
39 MeV/u for the Xe+Sn system. In these works, the nature of the transition 
could not be determined. 

We tested such a behaviour for the semi-peripheral collisions of interest in the 
present paper. In figures 23 and 24 the /(-2a) functions (with m = Z rnax ) are 
plotted for various Etrans bins for the 80 MeV/u Xe + Sn system. The plots 
correspond to A values of 0.5 and 1 respectively. No scaling is observed. How- 
ever, scaling is recovered if one separates the events according to the Z asym 
ranges corresponding to defined events: residues (Z asym larger than 0.7) or mul- 
tifragmentation (Z asym smaller than 0.3). The results are shown in figures 25 
and 26. Residue events are numerous only for Etrans values smaller than 200 
MeV. Only the corresponding bins have been retained in figure 25. Conversely, 
multifragmentation becomes dominant for Etrans larger than 200 MeV. Only 
the corresponding bins have been retained in figure 26. Similar results have 
been obtained for the Au+Au system with fission events reconstructed as 
explained in section 3 [49]. 

The fact that the scaling is clearly observed for the two classes of events sug- 
gests that these latter can be understood as an ordered (resp. disordered) 
phase. Conversely, the absence of scaling for the global distribution is consis- 
tent with the interpretation in terms of phase coexistence. Indeed, the absence 
of bimodality in the distribution of Z max may be due to a too strong corre- 
lation between the sorting variable and the deposited energy, which is also 
an order parameter: if a constraint is applied on an order parameter in the 
coexistence zone, its bimodal character is naturally suppressed. As the system 
is still in the coexistence zone, a pure phase scaling may be expected to fail 
as it is suggested in recent theoretical calculations [15]. 



7.2 Relation between bimodality and negative heat capacity 

Another signal of phase transition is negative heat capacity. This phenomenon 
is univocally associated to first order phase transitions with a finite latent heat 
in isolated systems with a size comparable to the range of the forces governing 
its equation of state [36]. A negative heat capacity occurs whenever the entropy 
presents an anomalous curvature as a function of the available energy. It has 
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been shown in reference [17] that heat capacity may be deduced from the 
measurement of the fluctuations in the sharing of the available energy between 
independent degrees of freedom. Since the heat capacity can be negative only 
in the microcanonical ensemble, from an experimental point of view events 
with a defined excitation energy have to be selected. For such a sorting, the 
total heat capacity may be approximately written: 

C«77^ (4) 



where C\ is the partial heat capacity for a subset of degrees of freedom (for 
instance momentum space), T is the temperature and a 2 is the variance as- 
sociated with the sharing of the total available energy between the degrees of 
freedom subset and the other ones. The total heat capacity becomes negative 
when these fluctuations become larger that the fluctuations associated to the 
canonical ensemble. 

This formalism has been used to study central Xe + Sn collisions from 32 to 
50 MeV/u [2], and peripheral Au + Au collisions at 35 MeV/u [1]. In both 
cases, fluctuations exceeding the canonical expectation were found. For the 
Xe + Sn central collisions, the signal is observed at 32 and 39 MeV/u but 
not at larger bombarding energies for which the system would have left the 
coexistence region. 

We have tried to correlate this signal with the bimodality one which is dis- 
cussed in this paper. For this purpose, events must be sorted according to the 
excitation energy of the QP, obtained by calorimetry on the QP source as ex- 
plained in section 6.1. We have retained only events with a total reconstructed 
source charge equal to the projectile charge within 10%. 

The negative heat capacity study has been performed in studying the sharing 
of the excitation energy between the kinetic and the potential energy parts 
at the freeze-out configuration. Rewriting equation 4 in normalizing the heat 
capacity to the source size gives: 

C ~ \ (5) 

n <A sou ><r 2 k v > 



o\ is the variance of the kinetic energy distribution; < A sou > is the average 
mass number of the source; Ck is the partial kinetic energy heat capacity per 
nucleon: 

C k = 6 

dT y J 



23 



The potential energy E p has been obtained from the mass balance corrected 
for the Coulomb part in assuming that all the LCP are emitted from the IMF 
in a second step after the freeze-out: 

Mjm F 

E P= m i~ m ^u) + E cou i (7) 

i=l 



M IMF is the IMF multiplicity (Z > 3); the mj are the IMF primary masses; 
m sou is the reconstructed QP mass; the Coulomb energy balance at freeze-out 
E cou i has been calculated from the Wigner-Seitz formula[50]. 

The temperature of relation 5 has been calculated from the total kinetic energy 
by using for the level density parameter aj = 4*-: 



MjMp o 

<E k >=< a t >T 2 + <-(M -1)>T 

i=l 



(8) 
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Fig. 28. Similar to figure 27 but in se- 
lecting the most compact QP events: 
see text. 
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Fig. 27. Evolution with the available 
excitation energy in MeV/u of the 
two terms of the denominator of rela- 
tion 5 (see text): partial heat capac- 
ity (open stars), and fluctuation term 
(black points). System Au + Au at 80 
MeV/u. 

From equation 5 it turns out that the heat capacity will be negative if the sec- 
ond term of the denominator (fluctuation term) exceeds the first one (partial 
heat capacity). These two terms are plotted in figure 27. Fluctuations are far 
below the canonical expectation whatever the excitation energy is. This result 
is quite general and has been observed for quasi-projectiles from the Xe + Sn 
and Au + Au systems at any bombarding energy. 
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At this point, one may recall that the negative heat capacity signal is much 
more fragile than the bimodality one. First, the evaluation of fluctuations 
is very sensitive to the details of the calorimetric hypotheses and side feed- 
ing corrections [51]. Moreover, any energy fraction which has not been shared 
among all the available degrees of freedom but is still stored in some selected 
ones will induce a decrease of the fluctuations we are looking at. In other 
words, negative heat capacity is not observed if dynamical effects freeze a 
sizeable part of the available energy in a few defined degrees of freedom remi- 
niscent of the initial direction of the beam. In order to decrease the influence 
of dynamical effects, we have tried to select QP events in which the memory 
of the entrance channel is minimized. We did such a selection in section 5.3 
in cutting the cos(6) distribution of figure 13). This cut does not modify sig- 
nificantly the signal. To ensure a stronger reduction of dynamical effects, we 
selected compact events, i.e. events for which the QP decay is not elongated 
in the momentum space. Compactness has been estimated from the relative 
velocity between the largest QP fragment and the mean velocity of the other 
QP fragments: 

Vrel — Vz max MjMF _! 



A larger value v rd is expected if intermediate velocity IMF contribute to the 
QP. We have retained only the events for which v re i was smaller than 3 cm/ns 
which is a value expected from Coulomb repulsion. The corresponding results 
are shown in figure 28. Only with this selection, abnormal fluctuations appear 
for excitation energies between 3 and 5 — 6 MeV/u. This is a general feature for 
the Au + Au system as shown in figure 29 while, at variance, the fluctuations 
remain below their canonical expectation for any excitation energy for the Xe 
+ Sn system (see figure 30). This difference between the two systems is not 
fully understood. 




Fig. 29. Systems Au + Au at 60 (left), 80 (middle) and 100 (right) MeV/u. Evo- 
lution with the available excitation energy of the two terms of the denominator of 
relation 5: partial heat capacity (open stars), and fluctuation term (black points). 
Non compact events (see text) have been rejected. 
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Fig. 30. Systems Xe + Sn at 65 (left), 80 (middle) and 100 (right) MeV/u. Evo- 
lution with the available excitation energy of the two terms of the denominator of 
relation 5: partial heat capacity (open stars), and fluctuation term (black points). 
Non compact events (see text) have been rejected. 
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Fig. 31. Similar to figure 28 for the Au + Au 80 MeV/u system but the events have 
been sorted according to Etrans. Each symbol corresponds to a selected Etrans 
bin of 0.25 MeV/u width. The fluctuation peak is observed in the bins associated 
with the open points and to the crosses covering an Etrans /u range from 0.65 to 
1.15 MeV/u) in which the bimodality transition has been observed in figure 11. 



Let us focus now on the Au + Au system for which fluctuations are more im- 
portant. One may first remark that the excitation energy range corresponding 
to the fluctuation peak is quite independant of the bombarding energy. In 
figure 31 we have connected this signal with the bimodality one. Events have 
been sorted according to the Etrans variable used previously for bimodality 
studies and only the most probable excitation energy bins have been retained 
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for each Etrans. Each symbol of figure 31 corresponds to an Etrans bin. The 
fluctuation peak corresponds to the third (open points) and fourth (crosses) 
bins, i.e. for Etrans/u values covering the range [0.65-1.15 MeV/u] in which 
the bimodality transition has been observed in figure 11. This coherence indi- 
cates that abnormal fluctuations and bimodality may sign the same physical 
phenomenon for the same set of events, but sorted in two different ways. 



8 Conclusion 

We studied two symmetrical systems Xe + Sn and Au + Au from 60 to 
100 MeV/u. We focused on semiperipheral collisions for which it is possible 
to isolate a QP contribution. The sorting of the collisions has been achieved 
from the LCP transverse energy on the QT side and we observed QP properties 
thus avoiding as much as possible auto-correlation due to sorting. We selected 
events for which most of the available charge on the QP side has been detected 
and identified. The distribution of the charge of the heaviest fragment and the 
asymmetry between the two heaviest fragments produced in each event was 
analysed. We have observed a two-peaked distribution for all systems for a 
temperature range around 5 MeV. One peak corresponds to events dominated 
by a residue and the other one to multifragmentation. Fission events in the Au 
case have been treated together with the residue-like events. If this behaviour 
is strongly suggestive of a phase transition, it does not allow to conclude about 
the nature of the phenomenon. 

We observed a scaling on the total mass of the system since the transverse en- 
ergy for which bimodality is observed is proportional to this mass for a defined 
bombarding energy per nucleon. On the other hand, we observed a transition 
transverse energy roughly proportional to the beam energy. We interpreted 
this effect as due to pre-equilibrium effects which increase with the bombard- 
ing energy. This statement is supported by the fact that the bimodality signal 
is observed for smaller transverse energy values if one selects events which are 
closer to equilibrium. This interpretation is still open to discussion. Neverthe- 
less, it turns out that many features are in agreement with an interpretation 
in terms of a phase transition in which the major parameter is the deposited 
energy. As a matter of fact, it has been established that the deposited (excita- 
tion) energies at the transition are larger for multifragmentation than for the 
residue-like events. At the same time, it has been shown that the temperatures 
involved in both cases are similar. This second conclusion is however weakened 
by the fact that the available thermometers are different for the two sets of 
events. 

A very interesting feature is that the bimodality signal can be correlated with 
defined scaling properties in the distribution of the heaviest fragment size 
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and the observation of important conflgurational energy fluctuations which 
we have tentatively associated at least for the Au+Au system to negative 
heat capacity. 

Altogether this set of results is in agreement with the expectations from a 
phase transition belonging to the liquid-gas universality class, and is well re- 
produced by a dynamical model. Many points need to be clarified. If the 
convexity can be seen in the Z max -Z asym representation, this convexity is not 
very pronounced and cannot be observed in the Z max distribution. This leaves 
the question of the order of the transition open. The relevant order parameter 
and the associated nature of the transition is not clear either. In particular the 
dynamical code HIPSE suggests that angular momentum may play a major 
role in the transition. 

Up to now, bimodality has been mainly recognized in semiperipheral collisions. 
The fact that it is more difficult to observe it in central ones may reflect a 
lower involved angular momentum. But it can also be due to the fact that the 
central collision sorting is too close to a microcanonical selection (for which 
bimodality should not be observed) rather than a canonical one ensuring the 
energy fluctuations which are necessary to observe the coexisting phases. In 
order to better elucidate this question, we intend to progress in two directions. 
We will first concentrate on a better understanding of central collisions in 
mixing various bombarding energies to simulate a canonical sorting. Another 
kind of analysis will concern the understanding of the role of the distance to 
equilibrium which is reached for a set of events. We will address this question in 
comparing collisions with an entrance channel which is asymmetric in isospin. 
The distance to equilibrium will be inferred from the distance to equilibrium 
in the isospin degree of freedom. Such measurements have been performed 
with INDRA at GSI for the 129 Xe + 112 Sn et 124 Xe + 124 Sn at 100 MeV/u. 
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